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HIGHLIGHTS 


•  New  charging  strategy  is  proposed  based  on  integration  of  Taguchi  method  and  stage  of  charge  estimation. 

•  Optimal  charging  current  pattern  is  determined  using  multiple-response  optimisation  technique  in  Taguchi  method. 

•  Experiments  of  lithium  ion  batteries  verify  the  effectiveness  of  the  new  charging  strategy. 

•  New  charging  strategy  has  better  performances  than  conventional  constant  current  constant  voltage  charging  method. 
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In  this  paper,  a  new  charging  strategy  of  lithium-polymer  batteries  (LiPBs)  has  been  proposed  based  on 
the  integration  of  Taguchi  method  (TM)  and  state  of  charge  estimation.  The  TM  is  applied  to  search  an 
optimal  charging  current  pattern.  An  adaptive  switching  gain  sliding  mode  observer  (ASGSMO)  is 
adopted  to  estimate  the  SOC  which  controls  and  terminates  the  charging  process.  The  experimental 
results  demonstrate  that  the  proposed  charging  strategy  can  successfully  charge  the  same  types  of  LiPBs 
with  different  capacities  and  cycle  life.  The  proposed  charging  strategy  also  provides  much  shorter 
charging  time,  narrower  temperature  variation  and  slightly  higher  energy  efficiency  than  the  equivalent 
constant  current  constant  voltage  charging  method. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  use  of  portable  electronic  devices  has  been  significantly 
increased  in  recent  years.  The  traditional  rechargeable  batteries 
such  as  lead-acid,  nickel-cadmium  and  nickel-metal  hydride  bat¬ 
teries  used  in  these  devices  have  tended  to  be  replaced  with 
lithium-ion  batteries,  which  exhibit  several  appealing  features  of 
high  power  and  energy  densities,  broad  operating  temperature 
range,  long  life  cycles,  no  memory  effect  and  low  self-discharge  rate 
[1-3  .  The  merits  of  these  lithium-ion  batteries  also  meet  high- 
energy  demand  of  electric  vehicles  (EVs),  making  them  as  very 
promising  primary  power  sources  for  EVs  in  the  future  [4  .  For  the 
applications  of  lithium-ion  batteries,  a  well-designed  battery 
charger  plays  a  vital  role  for  sustaining  battery  performance  and 
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lifespan  [5]  and  the  key  to  the  charger  is  the  charging  strategy 
including  the  selection  of  charging  current  pattern,  the  control  and 
termination  of  charging  process. 

The  simple  and  common  charging  strategy  of  a  constant  current 
constant  voltage  (CCCV)  is  widely  used  in  commercial  chargers  for 
lithium-ion  batteries  [6  .  It  has  two  charging  stages.  In  the  first 
constant  current  (CC)  stage,  the  battery  is  charged  at  a  constant 
current  until  battery  terminal  voltage  reaches  a  pre-defined  value 
(e.g.  4.2  V).  In  the  second  constant  voltage  (CV)  stage,  a  voltage  of 
the  pre-defined  value  is  maintained  until  the  charging  current  is 
gradually  dropping  to  the  pre-defined  value  (e.g.  0.05  C),  where 
0.05  C  represents  the  charging  current  in  C  rate  which  is  defined  as 
the  ratio  of  charging  current  to  nominal  capacity  (C)  of  a  battery.  In 
this  study,  0.05  C  is  chosen  as  a  cut-off  current  to  stop  charging 
process  for  the  CCCV  charging  method. 

Among  existing  charging  methods  [7],  the  multi-stage  constant 
current  (MCC)  charging  method  is  considered  as  a  high-quality 
charging  pattern  with  the  advantages  of  long  cycle  life,  high 
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charge/discharge  energy  efficiency  and  short  charging  time  [8]. 
Many  attempts  were  proposed  to  optimise  charging  current  pattern 
for  MCC  charging  methods.  It  is  regarded  as  a  combinatorial  opti¬ 
misation  problem.  Such  a  problem  is  difficult  to  be  solved  by 
exhaustive  search  [9  .  On  one  hand,  global  optimisation  techniques, 
such  as  genetic  algorithm  (GA),  ant  colony  system  (ACS)  and  par¬ 
ticle  swarm  optimisation  (PSO),  were  used  to  search  for  optimal 
charging  current  pattern  [10-12  .  On  the  other  hand,  Taguchi 
method  (TM)  was  utilised  and  applied  an  orthogonal  array  (OA)  to 
design  charging  current  patterns  so  that  only  a  small  number  of 
experiments  were  required  in  order  to  search  optimal  charging 
current  pattern  [13,14  .  For  these  MCC  charging  methods,  the  pre¬ 
defined  voltage  (e.g.  4.2  V)  is  used  to  switch  each  of  the  charging 
stages  as  well  as  a  cut-off  voltage  to  stop  charging  process. 

In  this  paper,  an  optimal  charging  strategy  based  on  the  MCC 
charging  method  has  been  developed  with  the  integration  of  the 
TM  and  the  state  of  charge  (SOC)  estimation.  The  TM  is  applied  to 
search  optimal  charging  current  pattern.  The  estimated  SOC  is  used 
to  control  and  terminate  the  charging  process.  The  experimental 
results  of  LiPBs  have  been  used  to  verify  the  effectiveness  of  the 
new  charging  strategy.  The  remaining  parts  of  this  paper  are 
organised  as  follows:  Section  2  experimentally  explores  the  in¬ 
fluences  of  various  charging  current  levels  on  charging  time, 
charging  efficiency  and  temperature  variation  during  charging 
process  to  assist  in  selecting  charging  current  candidates  to  form 
charging  patterns;  Section  3  presents  a  new  charging  approach. 
Section  4  shows  the  experimental  results  to  validate  the  proposed 
charging  method  by  comparing  it  with  the  equivalent  CCCV 
charging  method.  Section  5  concludes  the  paper. 

2.  Experimental 

2.1.  Battery  testing  platform 

A  battery  testing  platform  has  been  constructed  to  conduct  the 
experiments.  It  consists  of  a  programmable  power  supply,  a  pro¬ 
grammable  electronic  load  and  a  safety  control  switch  box.  The 
LabVIEW  software  program  has  been  used  to  implement  charging 
strategies  by  setting  different  charge  currents  (or  current  profiles), 
sample  experimental  data  and  store  the  data  into  the  computer  for 
further  analysis.  The  details  of  this  platform  can  be  found  in  Ref. 
[15]. 

Turnigy  lithium-polymer  batteries  (LiPBs)  with  the  nominal 
capacities  of  5.8  Ah  and  5.0  Ah  are  used  in  the  experiments,  which 
correspond  to  the  dimensions  of  136  mm  x  50  mm  x  10  mm  and 
135  mm  x  50  mm  x  9  mm  and  the  weights  of  143  g  and  130  g, 
respectively.  These  batteries  comprise  of  LiMn204  cathode  and  an 
artificial  graphite  anode.  They  are  usually  laminated  with  the  solid 
electrolyte  and  the  separator  sheets,  and  enclosed  in  a  foil-type 
case.  They  have  the  nominal  voltage  of  3.7  V  with  the  maximum 
voltage  of  4.2  V  for  charge  and  the  minimum  cut-off  voltage  of  2.7  V 
for  discharge. 

2.2.  Testing  of  a  LiPB  using  CCCV  charging  method 

A  LiPB  with  the  nominal  capacity  of  5.8  Ah  (Bl)  is  selected  to 
conduct  experiments.  The  CCCV  charging  method  is  adopted  to 
investigate  the  influences  of  various  charging  current  levels  on 
charging  efficiency,  charging  time  and  temperature  variation  dur¬ 
ing  charging  process.  According  to  battery  manufacturer's  data¬ 
sheet,  it  is  recommended  that  the  charging  current  of  0.3  C  is  for 
normal  charge  and  the  maximum  charging  current  of  2.0  C  is  for 
fast  charge.  Thus,  five  charging  currents  of  2.0  C,  1.5  C,  1.0  C,  0.7  C 
and  0.3  C  in  the  CC  stage  of  the  CCCV  charging  method  have  been 
chosen  to  charge  the  battery  at  room  temperature,  where  initially 


the  battery  is  fully  discharged  (e.g.  SOC  =  0%).  In  the  CV  stage,  the 
preset  voltage  of  4.2  V  is  maintained  to  charge  the  battery  until  the 
charging  current  is  gradually  reduced  to  0.05  C,  which  means  that 
the  battery  is  fully  charged  (e.g.  SOC  =  100%).  Then,  this  fully 
charged  LiPB  is  discharged  at  the  current  of  1.0  C  until  the  battery  is 
fully  discharged  with  the  cut-off  voltage  of  2.7  V.  For  each  of  these 
five  charging  and  discharging  processes,  the  experimental  data  of 
the  currents,  voltages  and  temperatures  are  collected  for  analysis. 
The  corresponding  results  are  shown  in  Fig.  1  and  Table  1,  respec¬ 
tively.  Based  on  the  observation  of  these  experimental  results,  a  few 
points  can  be  made  to  help  in  proposing  a  new  charging  strategy 
which  can  achieve  high  energy  efficiency,  short  charging  time,  and 
low  temperature  variation  in  Section  3.1.1. 

•  The  charging  current  of  2.0  C  is  allowed  to  charge  the  battery  up 
to  almost  90%  of  nominal  capacity  in  the  existing  arrangement 
of  the  CCCV  charging  methods.  Since  2.0  C  is  the  maximum 
charging  current  for  the  battery,  if  it  is  used  to  charge  the  battery 
all  the  time,  the  cycle  life  of  the  battery  will  be  decreased 
significantly  because  lithium  plating  can  be  formed  at  high 
charging  current  particularly  when  the  battery  is  in  high  SOCs 
and  high  temperature  [16,17].  Thus,  high  charging  current 
should  be  avoided  in  high  SOCs  to  minimise  its  effect  on  battery 
cycle  life  and  the  proposed  charging  strategy  will  select  the 
charging  current  based  on  the  estimated  SOC. 

•  The  higher  the  charging  current  is,  the  shorter  the  charging 
time,  the  higher  the  temperature  variation  and  the  lower  the 
energy  efficiency  are.  The  highest  charging  current  of  2.0  C  re¬ 
sults  in  the  lowest  energy  efficiency,  shortest  charging  time  and 
largest  temperature  variation  while  the  lowest  charging  current 
of  0.3  C  results  in  the  highest  energy  efficiency,  longest  charging 
time  and  lowest  temperature  variation.  Thus,  the  combination 
of  different  levels  of  charging  currents  will  be  adopted  in  the 
proposed  charging  strategy  to  achieve  the  balance  of  energy 
efficiency,  charging  time  and  temperature  variation. 

•  The  ampere  hour  (Ah)  efficiencies  for  all  the  charging  currents 
are  almost  100%,  which  indicates  that  this  battery  can  deliver  all 
the  ampere  hours  at  1.0  C  no  matter  what  charging  currents 
varying  from  0.3  C  to  2.0  C  are  adopted.  Thus,  the  Ah  efficiency 


Fig.  1.  Experimental  results  for  different  charging  currents  in  CC  stage  of  CCCV 
charging  methods. 
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Table  1 

Experimental  results  of  different  charging  currents  in  CC  stage  of  CCCV  charging 
method. 


Charging 

current 

Charging  time 

Temperature 

variation 

Ah  efficiency 

Energy 

efficiency 

2.0  C 

2775  s 

5.58  °C 

100% 

96.14% 

1.5  C 

3271  s 

4.31  °C 

100% 

96.46% 

1.0  C 

4325  s 

3.70  °C 

99.95% 

96.84% 

0.7  C 

5691  s 

3.40  °C 

100% 

97.40% 

0.3  C 

12,078  s 

3.31  °C 

99.80% 

97.51% 

will  not  be  considered  in  the  TM  to  search  optimal  charging 
current  pattern. 

•  The  influences  of  charging  current  levels  on  energy  efficiency, 
charging  time,  temperature  variation  are  only  appreciable  when 
the  charging  current  differences  are  0.5  C  (e.g.  from  2.0  C,  1.5  C 
to  1.0  C)  and  0.3  C  (from  1.0  C,  0.7  C  to  0.3  C)  from  one  CCCV 
charging  method  to  another,  respectively.  Thus,  charging  cur¬ 
rent  candidates  are  not  necessary  to  be  chosen  as  many  as  those 
in  the  existing  MCC  charging  methods  with  the  current  differ¬ 
ence  of  0.05  C  12—14],  leading  to  significant  simplification  of 
the  searching  process  for  optimal  charging  current  pattern. 


2.3.  Testing  ofLiPBs  using  current  pulse  charge 

As  mentioned  earlier,  the  estimated  SOC  controls  and  termi¬ 
nates  the  charging  process  in  the  new  charging  strategy.  Thus,  it  is 
important  to  accurately  estimate  the  SOC  during  the  charging 
process.  In  this  study,  a  relationship  between  the  open  circuit 
voltage  (OCV)  and  the  SOC  is  used  to  estimate  the  SOC.  To  evaluate 
the  capacity  and  ageing  effect  on  such  a  relationship  for  the  same 
type  of  batteries,  the  pulse  current  charge  (PCC)  test  has  been 
conducted  for  four  LiPBs  with  the  capacities  of  5.8  Ah  (B1  and  B2) 
and  5.0  Ah  (B3  and  B4),  respectively,  where  B1  and  B2  have 
different  numbers  of  cycles  and  so  are  B3  and  B4.  The  procedure  of 
conducting  the  PCC  test  is  explained  as  follows.  Initially,  four  LiPBs 
are  fully  discharged  (SOC  =  0%)  and  left  in  the  open  circuit  state  for 
an  hour  so  that  the  OCV  is  able  to  reach  its  equilibrium  state  at 
room  temperature.  Then,  they  are  charged  by  a  sequence  of  1.0  C 
pulse  current  consisting  of  six  minutes  charge  equivalent  to  10%  of 
the  SOC  increment  and  one  hour  rest  until  the  charge  current 
reaches  the  cut-off  current  of  0.05  C  when  the  batteries  are 


SOC  |%| 


considered  to  be  fully  charged  (SOC  =  100%).  The  experimental 
results  of  these  four  batteries  are  illustrated  in  Fig.  2.  It  can  be  seen 
that  the  relationship  between  the  OCV  and  the  SOC  for  these  four 
batteries  are  almost  the  same,  which  shows  that  the  SOC  estimation 
can  be  used  to  control  the  charging  process  for  these  types  of  LiPBs 
regardless  of  their  capacities  and  cycle  life. 

3.  Proposed  charging  strategy 

3.1.  Taguchi-based  optimal  charging  pattern  searching  algorithm 

The  TM  was  initially  developed  by  Genichi  Taguchi  for  opti¬ 
mising  manufacturing  process  and  then  accommodated  into  the 
fields  such  as  mechanical  engineering,  power  electronics,  and  radio 
networks  [18-21  .  For  the  optimisation  problems  with  all  possible 
combinations  of  parameters,  the  TM  is  suitable  for  finding  the 
optimal  parameters  and  decreasing  the  number  of  necessary 
experiments. 

3.1.1.  Charging  current  pattern  formulation 

According  to  the  discussions  in  Section  2.2,  a  four-stage  constant 
current  (FSCC)  charging  method  is  proposed  in  this  paper  and 
conceptually  shown  in  Fig.  3.  The  total  charging  period  is  equally 
divided  into  four  stages,  where  each  stage  has  the  range  of  the  SOC 
equal  to  25%.  In  each  stage,  the  LiPB  is  charged  by  a  pre-set  current. 
During  charging  process,  the  SOC  is  estimated.  When  the  estimated 
SOC  reaches  the  predetermined  SOC,  the  charging  process  will  be 
shifted  to  next  stage  and  a  new  pre-set  charging  current  will  be 
applied  accordingly.  The  charging  process  will  continue  until  the 
estimated  SOC  reaches  100%  when  the  battery  is  fully  charged  and 
the  charging  process  terminates. 

In  the  proposed  FSCC  charging  method,  the  combination  of 
different  levels  of  charging  current  in  four  stages  will  lead  to 
different  charging  current  patterns.  Fable  2  shows  the  charging 
current  candidates  in  this  study  where  the  charging  current  in  each 
stage  has  three  different  levels  with  the  current  difference  of  0.1  C, 
the  charging  current  in  different  stages  has  the  current  difference  of 
0.5  C  and  the  charging  currents  at  the  present  stage  are  lower  than 
those  of  previous  stage.  From  Table  2,  the  number  of  a  full  factorial 
charging  pattern  requires  81  experiments.  In  each  experiment,  a 
specific  charging  current  pattern  charges  the  fully  discharged  bat¬ 
tery  to  determine  the  charging  time  and  the  temperature  variation 
over  the  charging  process  and  then  the  battery  is  discharged  at  a 
predefined  constant  current  of  1.0  C  to  evaluate  the  influence  of 
each  charging  current  pattern  on  energy  efficiency.  Obviously,  it  is 
not  wise  to  test  all  81  possible  charging  current  patterns  to  search 


Fig.  2.  Experimental  results  for  relationship  of  OCV  versus  SOC  for  LiPBs. 


Fig.  3.  Proposed  charging  current  pattern. 
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Table  2 

Charging  current  candidates  for  four  charging  stages. 


Stage  name 

High 

Medium 

Low 

Stage  Si 

2.0  C 

1.9  C 

1.8  C 

Stage  S2 

1.5  C 

1.4  C 

1.3  C 

Stage  S3 

1.0  C 

0.9  C 

0.8  C 

Stage  S4 

0.5  C 

0.4  C 

0.3  C 

the  optimal  one.  To  eliminate  the  lengthy  experimental  procedure 
while  taking  advantage  of  the  performance  of  a  full  factorial 
experimental  design,  the  TM  is  developed  to  search  in  the  param¬ 
eter  space  with  a  small  number  of  experiments  based  on  the 
orthogonal  arrays  (OAs)  technique.  According  to  the  most 
commonly  used  standard  OAs  for  experimental  design  [14],  the 
search  of  an  optimal  charging  current  pattern  can  be  performed 
among  only  9  patterns  represented  by  L9  (34),  where  3  and  4  are  the 
number  of  charging  current  levels  in  each  stage  and  the  number  of 
stages  in  each  charging  current  pattern  shown  in  able  2,  respec¬ 
tively.  These  9  charging  current  patterns  as  shown  in  fable  3 
correspond  to  9  experiments.  With  these  patterns  selected  on  the 
basis  of  the  TM  method,  the  interactions  between  the  current  levels 
are  quasi-uniformly  distributed  between  the  columns  of  Table  3 
and  are  nearly  independent  [22  .  In  Fable  3,  the  High,  Medium 
and  Low  represent  different  levels  of  charging  current  for  each 
stage,  respectively.  For  example,  the  charging  current  pattern  of 
High,  High,  High  and  High  in  experiment  no.  1  stands  for  the  FSCC 
charging  pattern  with  the  charging  current  sequences  of  2.0  C,  1.5  C, 
1.0  C,  0.5  C.  Once  each  experiment  is  completed,  the  experimental 
results  are  considered  as  the  responses yy  with  i  =  1,  2,  ...,  9  which 
is  the  number  of  the  experiment  and  j  =  1,  2,  3  which  is  corre¬ 
sponding  to  energy  efficiency,  charging  time  and  temperature 
variation,  respectively. 

The  experiments  of  9  charging  current  patterns  are  conducted 
on  5.8  Ah  LiPB  (Bl)  at  room  temperature.  During  the  charging 
process,  the  SOC  estimation  is  carried  out  (see  the  detailed  expla¬ 
nation  of  the  SOC  estimation  in  Section  3.2).  When  the  estimated 
SOC  reaches  to  25%,  50%  and  75%,  respectively,  the  charging  process 
will  switch  to  the  next  stage.  When  the  estimated  SOC  reaches  to 
100%,  the  charging  process  terminates.  Once  the  charging  process  is 
completed,  the  battery  is  discharged  at  the  constant  current  of  1.0  C 
until  the  battery  voltage  reaches  to  the  cut-off  voltage  of  2.7  V. 
During  charging/discharging  processes,  if  the  battery  temperature 
is  higher  than  45  °C,  the  process  also  terminates  to  protect  the 
battery  from  thermal  runaway.  The  energy  efficiency  over  the 
charging  process  is  calculated,  and  the  charging  time  and  temper¬ 
ature  variation  are  collected.  All  responses  from  the  experimental 
results  are  shown  in  Table  4.  For  example,  the  experiment  no.  9 
(i  =  9)  has  three  responses:  energy  efficiency  {j  =  1)  of  95.96%, 
charging  time  (j  =  2)  of  3618  s  and  temperature  variation  (j  =  3)  of 
1.78  °C. 


Table  4 

Responses  (experimental  results)  for  energy  efficiency,  charging  time  and  temper¬ 
ature  variation. 


Exp.  No.  (i) 

Responses  (y,j) 

Energy  eff.  (%) 
for  j  =  1 

Charging  time  (s) 
for  j  =  2 

Temp.  var.  (°C) 
for  j  =  3 

1 

95.53 

3313 

4.97 

2 

95.12 

3939 

3.60 

3 

95.30 

4826 

2.99 

4 

95.11 

4571 

4.22 

5 

95.15 

3599 

2.97 

6 

95.11 

3932 

2.73 

7 

95.99 

4070 

3.95 

8 

95.68 

4602 

2.17 

9 

95.96 

3618 

1.78 

3.1.2.  Searching  optimal  charging  current  pattern 

The  purpose  of  the  new  charging  strategy  is  to  improve  charging 
efficiency,  reduce  charging  time  and  narrow  down  temperature 
variation  in  the  comparison  of  CCCV  charging  methods.  To  achieve 
this  by  using  the  TM,  the  experimental  results  shown  in  able  4  are 
divided  into  two  types:  the  larger-the-better  (LTB)  for  energy  effi¬ 
ciency  and  the  smaller-the-better  (STB)  for  both  charging  time  and 
temperature  variation.  The  TM  employs  signal  to  noise  (S/N)  ratio 
as  a  quality  measurement.  For  each  type  of  response  (j  =  1,  2  or  3), 
the  optimal  factor  level  is  the  factor  level  which  maximises  S/N 
ratio,  (3ij. 


For  LTB  type  response 


For  STB  type  response 


where  n  is  the  number  of  replicates  for  yp,  where  k  =  1,2,3,  ...,  n. 

In  order  to  solve  multi-responses  problem  in  the  TM,  the 
following  steps  are  utilised  [23]: 


Step  1:  The  S/N  ratios  (ft/)  for  energy  efficiency  (LTB  type), 
charging  time  and  temperature  variation  (STB  type)  are  calcu¬ 
lated  from  the  responses  (experimental  results)  from  Table  4 
using  Eq.  (1)  and  shown  in  Table  5. 

Step  2:  The  average  S/N  ratios,  ft  m  n,  for  three  types  of  responses 
(1  stands  for  energy  efficiency,  2  for  charging  time,  3  for  tem¬ 
perature  variation)  with  three  levels  (n  =  charging  current 
levels:  High,  Medium,  Low)  of  four  factors  (m  =  charging  stages: 
SI,  S2,  S3  and  S4)  are  estimated  from  responses  (ft/)  in  able  5. 


Table  3 

All  9  charging  current  patterns  based  on  L9  (34)  orthogonal  array. 


Exp.  no.  Stage  1 

Stage  2 

Stage  3 

Stage  4 

Responses 

O' =L2,  3) 

1 

High 

High 

High 

High 

yy 

2 

High 

Medium 

Medium 

Medium 

V2j 

3 

High 

Low 

Low 

Low 

V3 j 

4 

Medium 

High 

Medium 

Low 

Y4 j 

5 

Medium 

Medium 

Low 

High 

ys. j 

6 

Medium 

Low 

High 

Medium 

y  6j 

7 

Low 

High 

Low 

Medium 

V7 J 

8 

Low 

Medium 

High 

Low 

ys i 

9 

Low 

Low 

Medium 

High 

y  9/ 

Table  5 

S/N  ratios  for  energy  efficiency,  charging  time  and  temperature  variation. 


Exp.  No.  (i) 

Responses  (ft,-) 

Energy  eff.  (dB) 
for  j  =  1 

Charging  time  (dB) 
for  j  =  2 

Temp.  var.  (dB) 
for  j  =  3 

1 

39.603 

-70.404 

-13.925 

2 

39.565 

-71.908 

-11.118 

3 

39.582 

-73.672 

-9.516 

4 

39.565 

-73.200 

-12.505 

5 

39.568 

-71.124 

-9.448 

6 

39.565 

-71.892 

-8.732 

7 

39.645 

-72.192 

-11.926 

8 

39.616 

-73.259 

-6.740 

9 

39.642 

-71.169 

-5.023 
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To  illustrate  the  calculation  of  ftmn,  assume  that  ft^High  rep¬ 
resents  the  average  S/N  ratio  for  Si^ighl  it  is  an  average  of  the  S/N 
ratios  for  the  experiment  no.  1,  2,  and  3  for  energy  efficiency, 

namely  ftSl  ,High  —  3  (fti  +  fei  +  fti)  —  39.583  dB,  and 
ft ,sT ,Medium(39.566  dB)  and  ft,Sl,Low(39.634  dB)  are  the  average 
of  the  S/N  ratios  of  experiment  no.  4,  5,  6  and  experiment  no.  7, 
8,  9,  respectively.  In  the  same  wayftSi  High>  ft, .s, , Medium  and 
ft, St, Low  are  calculated  for  j  =  2  and  3,  and  the  average  S/N  ratios 
for  the  stages  2-4  can  also  be  calculated  from  each  response  ft)-. 
As  a  result,  all  the  average  S/N  ratios  are  shown  in  Table  6.  The 
TM  in  Refs.  [24,25  is  utilised  to  select  the  main  effects  for  each 
charging  stage.  According  to  Table  6,  the  optimal  factor  levels 
with  the  number  highlighted  in  bold  are  SUow,  S2,High,  $3, Low* 
^4,High;  for  energy  efficiency,  SiHigh*  SftHigh*  S3, High*  S^High  for 
charging  time  and  SUow*  $2, Low*  S3, Medium*  S4iHigh  for  temperature 
variation.  The  optimal  charging  pattern  is  the  balance  of  these 
main  effects  and  decided  by  using  average  level  weight  in  next 
step. 

Step  3:  Let  Wj>m>n  be  the  weight  of  level  n  for  factor  m  from 
response  ftmn,  which  is  calculated  as 


wj,m,n 


ft>,n 

maxnft-m,„ 

< 

maxnftm  n 
,  ft>,n 


For  LTB  type  response 
For  STB  type  response 


According  to  Eq.  ,  s^High >wi,Si, Medium  and  wi,Si,low  of  3  levels 
of  SI  are  calculated  for  energy  efficiency  (j  =  1 )  as  0.9987(=39.583 / 
39.634),  0.9983(^39.566/39.634),  1.0000(=39.634/39.634),  and 
w/,s1,High.w/,s1 , Medium *^j,Sl, Low  can  also  be  calculated  for  charging  time 
0  =  2)  and  temperature  variation  0  =  3),  respectively.  Similarly,  the 
calculations  of  all  wJ>m>n  can  be  made  for  stages  2-4,  they  are  all 

shown  in  Table  7.  Let  Wm  n  be  the  average  of  Wjim>n  over  all  responses 
at  the  same  level  (n)  and  factor  (m),  then  WSl  High,  WSl  Medium, 

WSl  low  are  calculated  from  three  responses  in  the  SI  as  0.8947, 
0.9231,  0.9990,  respectively.  Since  the  maximum  level  weight  for  SI 
corresponds  to  Si>Low.  the  optimal  level  of  SI  is  decided  as  Low. 
Similarly,  Wm,n  values  are  calculated  for  S2,  S3  and  S4,  where  their 
optimal  values  are  Low,  Medium,  and  High,  respectively.  Therefore, 
the  optimal  charging  current  pattern  is  the  sequence  of  1.8  C,  1.3  C, 
0.9  C  and  0.5  C. 


3.2.  Online  SOC  estimation 

Online  SOC  estimation  is  used  to  control  the  change  of  the 
charging  stage  and  terminate  charging  process  without 


Table  6 

Average  S/N  ratios  for  energy  efficiency,  charging  time  and  temperature  variation. 


Responses  (dB) 

Factor  (m) 

Stage  Si 

Stage  S2 

Stage  S3 

Stage  S4 

Level (n) 

Energy  eff. 

High 

39.583 

39.604 

39.595 

39.604 

Medium 

39.566 

39.583 

39.591 

39.591 

Low 

39.634 

39.596 

39.598 

39.588 

Charging  time 

High 

-71.995 

-71.932 

-71.852 

-70.899 

Medium 

-72.072 

-72.097 

-72.092 

-71.997 

Low 

-72.207 

-72.245 

-72.329 

-73.377 

Temp.  var. 

High 

-11.520 

-12.785 

-9.799 

-9.465 

Medium 

-10.228 

-9.102 

-9.549 

-10.592 

Low 

-7.896 

-7.757 

-10.297 

-9.587 

The  bold  number  indicates  the  biggest  number  in  each  stage,  the  level  corre¬ 
sponding  to  that  number  has  the  largest  effects. 


Table  7 

Level  weights  for  energy  efficiency,  charging  time  and  temperature  variation. 


Response  (dB) 

Factor  (m) 

Stage  Si 

Stage  S2 

Stage  S3 

Stage  S4 

Level (n) 

Energy  eff. 

High 

0.9987 

1.0000 

0.9999 

1.0000 

Medium 

0.9983 

0.9995 

0.9998 

0.9997 

Low 

1.0000 

0.9998 

1.0000 

0.9996 

Charging  time 

High 

1.0000 

1.0000 

1.0000 

1.0000 

Medium 

0.9989 

0.9977 

0.9967 

0.9847 

Low 

0.9971 

0.9957 

0.9934 

0.9662 

Temp.  var. 

High 

0.6854 

0.6067 

0.9745 

1.0000 

Medium 

0.7720 

0.8522 

1.0000 

0.8936 

Low 

1.0000 

1.0000 

0.9274 

0.9873 

Level  weight 

High 

0.8947 

0.8689 

0.9915 

1.0000 

Medium 

0.9231 

0.9498 

0.9988 

0.9593 

Low 

0.9990 

0.9985 

0.9736 

0.9844 

The  bold  number  indicates  the  selection  of  level  of  charging  current  for  each  stage. 


overcharging  LiPBs.  An  adaptive  switching  gain  sliding  mode 
observer  (ASGSMO)  based  on  battery  equivalent  circuit  model 
(BECM)  is  adopted  in  this  paper  for  the  SOC  estimation  15,26]. 

3.2.1.  Adaptive  switching  gain  sliding  mode  observer  for  SOC 
estimation 

A  variety  of  battery  models  has  been  developed  to  capture  dy¬ 
namic  characteristics  of  LiPBs  for  achieving  accurate  estimation  of 
battery  states  [27-29].  A  BECM  consisting  of  lumped  circuit  com¬ 
ponents  such  as  capacitors,  resistors,  voltage  sources  is  utilised  in 
this  paper.  The  schematic  diagram  of  the  BECM  is  shown  in  Fig.  4, 
where  the  variations  of  the  circuit  parameters  are  included  to 
represent  the  modelling  errors. 

In  the  BECM,  a  capacitor  Cn  as  the  nominal  capacity  represents 
the  total  stored  energy  in  the  battery.  The  current  source,  I  denotes 
the  charge  or  discharge  current  of  the  LiPB  and  Vt  expresses  the 
corresponding  battery  terminal  voltage.  The  voltage  across  the  Cn  as 
the  OCV,  Voc  varies  in  the  range  of  the  SOC,  Z  from  0%  to  100%  and  it 
represents  the  SOC  quantitatively.  An  ohmic  resistor  Rin  charac¬ 
terises  the  charging  energy  losses  of  battery.  Two  parallel  resistive 
and  capacitive  networks  comprise  of  the  electrochemical  polar¬ 
isation  resistance  Rpe,  capacitance  Cpe,  the  concentration  polar¬ 
isation  resistance  Rpc  and  capacitance  Cpc,  which  reflect  the  short¬ 
term  and  long-term  transient  responses  due  to  the  relaxation  ef¬ 
fect  of  battery,  respectively.  This  relaxation  effect  represents  the 
slow  convergence  of  battery  terminal  voltage  to  the  OCV  at  its 
equilibrium  state.  The  voltage-controlled  voltage  source,  Voc  (Z)  is 
used  to  bridge  the  nonlinear  relationship  between  the  SOC  and  the 
OCV  as  shown  in  Fig.  2.  The  self-discharge  resistance  is  ignored  in 
this  model  as  the  LiPB  has  extremely  low  self-discharge  rate.  The 
symbols  AV0O  ARin,  ARpe,  ACpe,  ARpc  and  ACpc  reflect  variations  of 
the  circuit  parameters  caused  by  the  adoption  of  the  constant 
values  of  the  BECM  parameters  and  the  piecewisely  linearised 
relationship  of  OCV-SOC  curves. 


Rpe  &ARpe  Rpc  &  ARpc 


Fig.  4.  Schematic  of  battery  equivalent  circuit  model. 
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The  SOC  indicates  the  ratio  of  charged  capacity  to  the  nominal 
capacity,  and  it  is  expressed  as. 


x(t)  =  Ax(t)  +  Bu(t )  +  A f(x,  u,  t) 


(ii) 


t 

Z(t)=Z(0)  +  (l/C„)  J  vl(r)dr  (3) 

0 

where  Z(0)  denotes  the  initial  SOC,  Cn  is  the  battery  nominal  ca¬ 
pacity  in  Ah,  J(r)  is  the  instantaneous  charge  current  (assumed  to  be 
positive),  and  r\  is  Ah  efficiency,  which  is  close  to  100%  for  LiPBs  in 
the  wide  range  of  current  and  temperature  as  shown  in  able  1. 

According  to  Kirchhoff  voltage  law  in  circuit  theory,  the  battery 
terminal  voltage  Vt  can  be  expressed  as. 


y(t)  =  Cx(t)  (12) 

where 


~-a-i 

a-i 

-a3 

-a4“ 

~~b\~ 

A  = 

ci2 

-a2 

02 

02 

,B  = 

0 

0 

0 

-a4 

0 

b2 

0 

0 

0 

-03- 

-  ^3  _ 

and  x(t) 

=  [Vt 

Z 

Vpe 

e- - 

1 — 1 

u 

cx 

Vt  —  Voc{Z)  +  Vpe  +  VpC  +  IR[n  +  Auncerta}n  (4) 

The  time  derivatives  of  the  SOC  and  polarisation  voltages  yield. 

Z=(//C„)  +  A/2  (5) 

\/pe  =  —  Vpe/(^peCpe)  +  f/Cpe  +  A/3  (6) 

V  pc  —  —  VpC/(i^pcCpc)  +  //Cpc  +  A/4  (7) 


where  Vpe  and  Vpc  denote  the  electrochemical  and  concentration 
polarisation  voltages  across  Cpe  and  Cpc,  respectively.  The  terms 
^uncertain.  A/2,  A/3  and  A/4  represent  model  uncertainties  and  time- 
varying  internal  and  external  disturbances. 

Despite  the  nonlinearity  of  OCV-SOC  curves  as  shown  in  Fig.  2, 
there  exists  a  piecewise  linear  relationship  between  the  OCV  and 
the  SOC  in  a  certain  range  of  the  SOC.  Thus,  the  OCV  can  be 
expressed  as  a  function  of  the  SOC  by  using  piecewise  linearisation 
method. 

Voc(Z)  =  kZ  +  v  (8) 

where  the  linearised  coefficients  of  k  and  v  are  varying  in  different 
SOC  ranges. 

Due  to  the  fast  sampling  rate  and  constant  charging  current 
rates  in  each  stage,  the  time  derivative  of  charging  current  can  be 
negligible,  namely  dl/dt  =  0,  and  thereby  the  time  derivative  of 
terminal  voltage  Vt  in  Eq.  (4)  with  the  substitutions  of  Eqs.  (5)— (7) 
gives. 


The  unknown  function  A/x,u,t)  represents  the  matched  model 
uncertainties  and  it  can  be  assumed  that. 

A =  r$(x,t)  (13) 

where  r  is  the  model  uncertainties  input  matrix  and  the  function 
?(x,t)  is  unknown  but  bounded 

|f(x,t)|  <  V'  VxeR4,  t>  0  (14) 

where  the  bound  /  is  a  positive  constant  and  can  be  determined  by 
the  largest  modelling  errors  between  the  experimental  data  and 
the  BECM  of  the  LiPBs.  Substituting  Eq.  (13)  into  Eq.  (11)  gives 

x(t)  =  Ax(t)  +  Bu(t)  +  r?(x,  t)  (15) 

According  to  [25],  the  ASGSMO  for  the  SOC  estimation  of  LiPBs  is 
derived  as  follows. 

x(t)  =  Ax(t)  +  Buft )  +  K(y(t)  -  Cx(t))  +  v  (16) 

y(t)  =  Cx(t)  (17) 

where  the  feedback  gain  matrix,  I<  and  the  adaptive  switching  gains 
function,  v  can  be  designed  so  that  the  stability  and  robustness  of 
the  ASGSMO  is  ensured. 

?(ey(t))rsgn(ey(f))  if  ey(t)*0  (18) 

0  otherwise 

The  discontinuous  switching  term  is  defined  as. 


V t  —  k(1  / Cn)  —  Vpc!  (RpeCpe)  +  // Cpe  —  Vpc/ (/^pcCpc)  +  // (-pc  +  A f\ 

(9) 

Solving  /  in  Eq.  (4)  and  substituting  it  into  Eq.  (5)  as  well  as  re¬ 
arranging  Eqs.  (6) — (8)  yield  the  state-space  equations  of  the  BECM 
as  follows. 


Vt  =  -at  Vt  +  at  Voc  (z)  +  a3Vpe  +  a4Vpc  +  bU  +  Afi 
Z  =  d'Wt  -  O'jVoc  (z)  —  a 2  Vpe  -  Cl3  Vpc  +  Af2  (10) 

Vpe  =  — ^Vpe  +  +  A/3 

Vpc  —  — ^3  VpC  -f-  +  A/| 

where  cq  =  l/(^pe6pe)  +  l/(^pc6pc)>  ^2  —  V(^inCn)*  ^3  —  V(^pc6pc)» 
a4  =  V(^pe6pe)»  =  K/Cn  +  R[n/(RpeCpe )  +  1/Cpe  +  R[n/{RpcCpc )  +  V 
Cpc,  ^2  —  VCpe  and  =  l/CpC. 

If  the  input  and  the  output  of  circuit  model  are  defined  as  u(t)  =  / 
and  y(t)  =  Vt,  respectively,  and  the  system  state  variables  are  chosen 
as  Vt,  Voc,  Vpe  and  Vpc,  then  Eq.  (10)  can  be  concisely  expressed  in 
matrices. 


sgn(ey(t)) 


+  1  ,  6y  ft)  >  0 

-1,  ey(t)<  0 


(19) 


where  ey(t)  =y(t) -y(t)  =  Ce(t)  and  the  adaptive  switching 
gain  d(Ceft))  is  updated  by  the  following  adaptive  law 


?(ey(t))  =  a\  ey(t) 


(20) 


where  a  is  a  positive  constant,  which  can  be  chosen  to  adjust  the 
speed  adaptation  for  switching  gains  function.  As  the  switching 
gains  are  updated  correspondingly  to  the  state  error  dynamics,  the 
ASGSMO  can  provide  robust  tracking  capability  against  modelling 
errors  while  significantly  reducing  the  magnitude  of  chattering  in 
the  SOC  estimation. 


3.2.2.  Battery  equivalent  circuit  model  parameters  identification 
The  parameters  of  the  BECM  shown  in  Fig.  4  are  extracted  from 
the  least  square  estimation  method  by  executing  a  PCC  profile  as  a 
battery  characterisation  test  at  room  temperature.  As  an  example, 
the  PCC  profile  is  performed  to  identify  the  BECM  parameters  of 
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5.8  Ah  LiPB  (B1 ).  The  corresponding  current  and  voltage  response  is 
shown  in  Fig.  5.  It  can  be  seen  that  there  are  totally  ten  sets  of 
transient  responses  as  the  battery  terminal  voltage  converged  to  its 
associated  OCV  after  one  hour  relaxation. 

From  Eq.  (4),  the  voltage  difference  between  terminal  voltage 
and  the  OCV  can  be  defined  as  Vd  =  (Vt-Voc)  such  that. 

Vd  =  VPe  +  Vpc  +  /Kin  (21) 

Based  on  Eqs.  (6),  (7)  and  (21),  the  system  transfer  function  of 
the  BECM  can  be  derived  as  follows. 


Eq.  (24)  can  be  further  concisely  expressed  in  following  form. 
U=0'0  (25) 

where  U  =  [ Vt(k  +  2)  Vt(k  +  3)  •••  Vt(n)]r, 

~-Vt(k  + 1)  -Vt(k)  I(k  +  2)  /(fc+1)  /(/<)  1" 

0=  —Vt(k  +  2)  Vt(k  +  1)  I(k  +  3)  I(k  +  2)  J(k+1)  1 

_  —Vt{n  -  2)  -Vt(n-1)  /(n)  /(n  —  1)  J(n  -  2)  1 


Vd(5) 


1? 


pe 


1  +  5i^peCpe 


+ 


1  +  Sl^pc 


(22) 


It  can  be  transformed  into  the  difference  equation  form. 


and  ©  =  [0\  d2  63  64  65  66}t  where  []T  denotes  the  trans¬ 

position  of  a  matrix.  Therefore,  the  least  square  method  can  be 
applied  to  estimate  the  parameters  based  on  the  following 
calculation 


Vt(k  +  2)  +  0\Vt(k  +  1)  +  02Vt{k) 

=  03I(k  +  2)  +  64 /(/<  +  1 )  +  d5I(k)  +  Oq 


(23) 


®=(&T-&)  -0TU 


(26) 


where  k  denotes  the  time  step  and  0\  =  l/(ftpeCpe)  +  l/(/?pCCpc), 
O2  —  l/(^pe6pe^pc6pc)>  ^3  —  ^in*  ^4  —  (^in  +  ^pe)/(^pe6pe)  +  (^in  +  ^pc)/ 
(i?pCCpC),  0$  =  ( R\n  +  i^pe  +  ^pc)/(^pe^pc6peCpc)  and  =  Voc,  and  the 
BECM  parameters  will  be  estimated  based  on  coefficients  of  Eq. 

(23). 

There  exist  n  sets  of  sampled  data,  as  shown  in  the  following 
matrix  form. 


where  ()_1  denotes  the  inverse  of  a  matrix.  Using  Eq.  (26),  the 
parameters  of  the  BECM  for  the  5.8  Ah  LiPB  (B1 )  are  Rm  =  3.15  mQ, 
ftpe  =  1.651  mfi,  RpC  =  2.35  mQ,  Cpe  =  1.7896  1<F  and 
Cpc  =  5.7368  1<F. 

To  verify  the  accuracy  of  the  BECM  with  the  identified  param¬ 
eters,  the  PCC  profile  is  reloaded  into  the  BECM.  Fig.  6  illustrates  the 
comparison  of  the  terminal  voltages  obtained  from  the  experiment 
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Fig.  5.  Pulse  current  charge  and  terminal  voltages  for  LiPB  Bl. 


Fig.  6.  Comparison  of  terminal  voltage  between  experimental  results  of  LiPB  Bl  and 
those  calculated  from  BECM  at  PCC  with  their  absolute  errors. 
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Time  [s] 

Fig.  7.  Comparison  of  terminal  voltage  between  experimental  results  of  LiPB  B1  and 
those  calculated  from  BECM  at  CCCV  with  their  absolute  errors. 

of  the  5.8  Ah  LiPB  (B1 )  and  those  calculated  from  the  BECM  as  well 
as  the  associated  modelling  errors.  It  shows  that  the  errors  are 
within  a  small  range  of  -0.025  V~  +  0.025  V.  Although  the  battery 
terminal  voltage  of  the  BECM  has  a  small  deviation  from  those  of  a 
real  LiPB  in  the  PCC  profile,  there  exist  considerable  modelling  er¬ 
rors  for  other  charging  current  profiles  when  the  predetermined 
constant  model  parameters  are  used.  For  example,  Fig.  7  shows  the 
comparison  of  the  terminal  voltage  obtained  from  the  experiment 
of  the  5.8  Ah  LiPB  (Bl)  and  those  calculated  from  the  BECM  for  the 
CCCV  charging  profile  with  1.0  C  charging  current.  Nevertheless,  the 
ASGSMO  for  the  SOC  estimation  can  overcome  the  inaccuracy  of  the 
BECM  with  the  constant  model  parameters  and  track  the  SOC 
accurately. 


SOC  directly  generated  by  the  Ah  counting  module  to  show  the 
accuracy  of  the  SOC  estimation.  The  feedback  gain  matrix,  I<  is 
obtained  by  using  the  LQR  method  with  Riccati  equation  as 
K=  [1.6872  2.7983  0.8271  0.3754 ]r  and  the  model  un- 
certainties  input  matrix,  T  is  defined  as  r  = 
[1.4516  -0.0326  0.0122  0.0204 ]r.  The  positive  constant,  a  is 
selected  as  0.4  to  satisfy  the  adaptation  speed  for  the  switching 
gains  function. 

The  first  charging  profile  is  the  PCC  and  its  corresponding  results 
of  the  SOC  estimation  from  charging  current  of  1.0  C  in  each  pulse 
charge  are  shown  in  Fig.  9.  It  can  be  seen  that  during  the  one  hour 
resting  phase  as  the  true  SOC  remains  same,  and  the  ASGSMO  is  able 
to  track  the  true  SOC  with  no  chattering  ripples,  also  its  SOC  esti¬ 
mation  errors  expressed  as  the  root  mean  square  errors  (RMSEs)  are 
limited  to  the  small  range  of  -4%  to  +2%.  This  is  due  to  the  fact  that 
the  switching  gains  of  the  ASGSMO  can  be  adaptively  tuned  to  the 
proper  low  levels  as  the  corresponding  estimation  errors  decrease. 

To  further  evaluate  the  ASGSMO  for  the  SOC  estimation,  another 
estimation  results  as  shown  in  Fig.  10  are  collected  from  the  CCCV 
profile  with  the  charging  current  of  1.0  C.  It  can  be  found  that  in  the 
presence  of  modelling  errors  caused  by  constant  model  parameters, 
the  ASGSMO  is  capable  of  tracking  true  SOC  with  the  RMSEs  in 
range  of  -3%  to  +3%.  The  estimated  SOC  is  rapidly  increasing  at  the 
CC  stage  and  it  becomes  slow  increment  at  the  CV  stage  as  the 
contant  battery  terminal  voltage  leads  to  the  quick  drop  in  charging 
current. 

4.  Results  and  discussions 

The  experiments  have  been  conducted  to  verify  the  effective¬ 
ness  of  the  new  charging  strategy,  which  has  the  charging  current 
sequence  of  1.8  C,  1.3  C,  0.9  C  and  0.5  C  and  is  used  to  charge  LiPB 
5.8  Ah  (Bl )  with  different  initial  SOCs  and  the  other  three  LiPBs  (B2, 
B3,  and  B4)  with  the  same  initial  SOC.  Furthermore,  the  new 
charging  strategy  compares  with  the  equivalent  CCCV  charging 
method.  The  followings  are  the  experimental  results  and  the  cor¬ 
responding  discussions. 

4.2.  Testing  results  on  LiPB  5.8  Ah  (Bl)  with  different  initial  SOCs 


3.2.3.  Verification  of  ASGSMO  for  SOC  estimation 

To  verify  the  effectiveness  and  accuracy  of  the  ASGSMO  for  the 
SOC  estimation,  two  charging  current  profiles  have  been  loaded  to 
the  5.8  Ah  LiPB  (Bl )  at  the  room  temperature  and  these  two  current 
profiles  are  simultaneously  applied  to  the  LiPB  integrated  with  Ah 
counting  module  and  the  ASGSMO  module  as  shown  in  Fig.  8.  The 
LiPB  terminal  voltage  is  sampled  and  fed  into  the  ASGSMO  to 
generate  the  output  tracking  error,  which  can  be  used  to  update  the 
corresponding  switching  gains  of  ASGSMO  for  the  compensation  of 
the  modelling  errors.  The  resultant  output  of  the  ASGSMO  module 
is  the  estimated  SOC,  which  is  concurrently  compared  to  the  true 


The  new  charging  strategy  has  been  used  to  charge  LiPB  5.8  Ah  (Bl ) 
at  different  initial  SOCs.  The  procedure  to  conduct  the  experiments  is 
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Fig.  8.  Verification  of  ASGSMO  for  SOC  estimation. 


Fig.  9.  ASGSMO  estimated  SOC  at  PCC  with  RMSE(%). 
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Fig.  10.  ASGSMO  estimated  SOC  at  CCCV  with  RMSE(%). 


explained  as  follows.  For  SOC  =  0%,  the  battery  is  fully  discharged  and 
rested  for  an  hour,  and  then  charged  by  the  new  charging  strategy.  For 
the  other  SOCs,  the  battery  is  fully  discharged  and  rested  for  an  hour, 
and  then  charged  by  the  constant  current  of  1.0  C  to  the  preset  SOCs, 
such  as  17%,  37%,  67%  and  84%  which  are  in  the  ranges  of  0%-25%, 
25%-50%,  50%-75%  and  75%-100%,  respectively,  and  then  charged 
by  the  new  charging  strategy.  Fig.  11  shows  the  charging  current 
patterns  and  the  corresponding  terminal  voltages.  It  indicates  that  the 
new  charging  strategy  can  successfully  charge  the  LiPB  at  any  initial 
SOCs  with  the  assistant  of  the  ASGSMO  for  the  SOC  estimation  which 
is  used  to  control  and  terminate  the  charging  process. 

4.2.  Testing  results  on  LiPBs  5.8  Ah  (B2),  5.0  Ah  (B3  and  B4)  at  the 
same  initial  SOC 

The  new  charging  method  has  also  been  used  to  charge  the 
other  three  LiPBs  5.8  Ah  (B2),  5.0  Ah  (B3  and  B4).  To  make 


±  15  n 
10 


G 

-  I 

u  5  1 

a  0 
U 


15 
c  10 

i  5 

a  0 


^  15 
c  10 
Z  5 

s_ 

a  0 


i.  15  n 


10  ■ 
5  ■ 


a> 

J  ° 


5 

4 


(a) 


Ho 


(b) 


(c) 


4.5 
4 

3.5 

4.5 
4 

3.5 

4.5 
4 

3.5 


&X 

£ 


o 

Ct 

£ 


SJ 

GJD 

£ 


CJ 

OX) 

Cl 


(d) 


40  60 

SOC  [%| 

(e) 


80 


100 


4.5  > 


01 

ox 

C5 

-w 

I 


Fig.  11.  Experimental  results  of  new  charging  strategy  with  different  SOCs  of  LiPB  B1 
(a)  SOC  =  0%,  (b)  SOC  =  17%,  (c)  SOC  =  37%,  (d)  SOC  =  67%,  (e)  SOC  =  84%. 
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Fig.  12.  Experimental  results  of  new  charging  strategy  at  SOC  =  0%  with  different 
LiPBs:  B2  (5.8  Ah),  B3  (5.0  Ah)  and  B4  (5.0  Ah). 


experimental  results  understand  easily,  the  initial  SOCs  for  these 
batteries  are  set  to  0%.  Fig.  12  shows  the  charging  current  pattern, 
terminal  voltage  and  temperature  variation,  their  voltage  and 
temperature  profiles  are  quite  consistent  for  these  three  batteries.  It 
shows  that  the  new  charging  strategy  can  apply  to  the  same  type  of 
LiPBs  with  different  capacities  and  cycle  life. 


Fig.  13.  Comparison  between  new  charging  strategy  and  CCCV  charging  method  on 
LiPB  Bl. 
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Table  8 

Experimental  comparisons  between  proposed  and  CCCV  charging  method  for 
different  LiPBs. 


LiPBs 

Charging 

algorithm 

Charged 

capacity 

(Ah) 

Discharged 

capacity 

(Ah) 

Energy 
eff.  (%) 

Charging 
time  (s) 

Temp, 
var.  (°C) 

Bl  5.8  Ah 

Proposed 

5.75 

5.66 

96.22 

3915 

1.73 

CCCV 

5.77 

5.76 

95.38 

4838 

3.4 

B2  5.8  Ah 

Proposed 

5.55 

5.49 

95.39 

3672 

2.05 

CCCV 

5.65 

5.64 

94.97 

4768 

3.27 

B3  5.0  Ah 

Proposed 

4.74 

4.64 

95.57 

3602 

2.61 

CCCV 

4.69 

4.65 

95.50 

4561 

3.18 

B4  5.0  Ah 

Proposed 

4.76 

4.66 

94.30 

3624 

1.37 

CCCV 

4.81 

4.80 

93.90 

5162 

3.21 

4.3.  Comparison  with  equivalent  CCCV  charging  method 

The  experimental  results  of  the  new  charging  strategy  are 
compared  with  those  of  the  equivalent  CCCV  charging  method  to 
illustrate  the  performance.  For  the  LiPBs  (Bl,  B2,  B3  and  B4),  the 
charging  currents  of  the  CC  stage  in  the  CCCV  charging  method  are 
chosen  as  the  average  charging  currents  of  the  new  charging 
strategy  correspondingly  and  calculated  form  Eq.  (27)  which  are 
0.910  C,  0.936  C,  0.944  C  and  0.942  C,  respectively. 

qq  _  l(^n  X  Tm) 

£m=l  Tm 

where  CC  is  the  average  charging  current,  Im  is  the  charging  current 
in  Tm  charging  duration  for  each  charging  stage  in  the  new  charging 
strategy. 

The  experimental  results  for  5.8  Ah  LiPB  (Bl)  are  shown  in 
Fig.  13.  The  charging  time  of  the  new  charging  strategy  is  shorter 
than  that  of  the  CCCV  charging  method  (80  min)  by  15  min,  which 
is  equivalent  to  22.5%.  The  temperature  variation  of  the  new 
charging  strategy  is  almost  half  of  that  of  the  CCCV  method.  To 
further  demonstrate  the  merits  of  the  new  charging  strategy,  the 
experiments  of  the  other  three  LiPBs  (B2,  B3  and  B4)  are  also 
conducted  and  the  corresponding  results  are  summarised  in 
Table  8.  It  can  be  seen  that  the  new  charging  strategy  achieves  the 
similar  better  results  to  those  obtained  from  Bl  than  the  corre¬ 
sponding  CCCV  charging  methods. 

4.4.  Application  of  proposed  charging  approach  to  other  types  of  Li- 
ion  batteries 

The  proposed  charging  strategy  has  been  developed  based  on 
the  observation  of  the  experimental  results  as  demonstrated  in 
Section  2.  If  the  other  types  of  Li-ion  batteries  have  the  similar 
experimental  results  from  the  CCCV  charging  tests  and  the  similar 
characteristics  of  the  OCV  versus  the  SOC  to  the  batteries  in  this 
study,  the  proposed  charging  strategy  should  be  able  to  apply  for 
these  types  of  Li-ion  batteries. 

5.  Conclusions 

The  new  charging  strategy  of  the  LiPBs  is  proposed  based  on  the 
integration  of  Taguchi  method  (TM)  and  the  SOC  estimation.  In  the 
TM,  the  principle  of  the  selection  of  charging  current  candidates, 
the  formation  of  charging  current  patterns  and  the  determination 


of  the  optimal  charging  current  pattern  is  elaborated.  In  the  SOC 
estimation,  the  procedure  to  design  an  adaptive  switching  gain 
sliding  mode  observer  (ASGSMO)  is  explained  in  detail.  The  effec¬ 
tiveness  of  the  new  charging  strategy  has  been  verified  by  the 
experimental  results  of  LiPBs  and  compared  with  the  equivalent 
CCCV  charging  method.  It  shows  that  the  new  charging  strategy  can 
successfully  charge  the  LiPB  with  different  initial  state  of  charges 
and  the  same  type  of  LiPBs  at  different  capacities  and  cycle  life. 
Furthermore,  it  has  the  better  performances  (much  shorter 
charging  time,  narrower  temperature  variation  and  slightly  higher 
energy  efficiency)  than  the  equivalent  CCCV  charging  method. 
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